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Abstract 

We consider gluodynamics in case when both color and magnetic charges are present. 
We discuss first short distance physics, where only the fundamental \Q m \ = 1 
monopoles introduced via the 't Hooft loop can be considered consistently. We show 
that at short distances the external monopoles interact as pure Abelian objects. 
This result can be reproduced by a Zwanziger-type Lagrangian with an Abelian 
dual gluon. We introduce also an effective dual gluodynamics which might be a 
valid approximation at distances where the monopoles \Q m \ = 2 can be considered 
as point-like as well. Assuming the monopole condensation we arrive at a model 
which is reminiscent in some respect of the Abelian Higgs model but, unlike the 
latter leaves space for the Casimir scaling. 



1 Introduction. 



Constructing gauge theories where the same gauge bosons interact with both magnetic 
end electric charges is a long standing problem. In case of electrodynamics, a particular 
solution is provided by the Zwanziger formalism |I] . In this formalism, one introduces for- 
mally two different gauge fields, A^, interacting with the electric and magnetic charges, 
respectively. To avoid the doubling of the degrees of freedom one imposes, however, the 
constraint that the field strength tensors build up on the potentials A^ and are dual 
to each other: 

F^A) = *F^ B ) = \^ Xp F Xp (B). (1) 

If one resolves the constraint and works in terms of a single potential, A^, then the string 
singularities naturally emerge. 

A direct evaluation of the radiative corrections within the Zwanziger formalism brings 
in a paradoxical conclusion that the Dirac quantization condition for the product of 
the electric and magnetic coupling, e ■ g = 2ir, gets violated (for a review and further 
references see H). The origin of this mismatch can be traced back to violation, within 
the perturbative approach of the Dirac veto which forbids direct interaction of particles 
with the Dirac strings. This turns out to be not true for virtual particles. There exist 
approaches, non-perturbative in nature, which respect the Dirac veto and then one can 
see, indeed, that the Dirac quantization condition survives the radiative corrections, for 
reviews see 0, |3|. At this point, a more detailed analysis of the Dirac strings is required 
so that one leaves, in a way, the ground of pure field theories. 

Our interest here is in constructing gluodynamics which would incorporate interaction 
both with color and magnetic charges. The history of this topic is also rich, see, e.g., ||, fj]. 
|5|, [| |7J and references therein. It is worth emphasizing from the very beginning that there 
are two different types of monopoles which can be considered in case of gluodynamics. 
Monopoles with the minimal magnetic charge, \Q m \ = 1 can be introduced via the 't Hooft 
loop ||. In the classical limit they correspond to classical solution with infinite action 
which are nothing else but the Dirac monopoles. These monopoles can be considered only 
as external sources to probe the vacuum state of the gluodynamics. The monopoles with 
\Q m \ — 2, on the other hand, do not exist as stable classical solutions ||. Nevertheless just 
these monopoles were observed and studied in great detail in the lattice simulations, for 



review and further references see, e.g., ||10|| . A crucial step to define \Q m \ = 2 monopoles 
is the partial gauge fixing, which leaves a remaining U(l) invariance JET]]. Then the 
monopoles are defined with respect to this U{\). 

The studies of \Q m \ = 2 monopoles strongly support the dual-superconductor mecha- 
nism of confinement |12| for the heavy quarks in the fundamental representation. Namely, 



the quark potential grows linearly at large distances: 

lim V q q(r) = or , (2) 



+ OO 



and the monopole contribution to V qq -(r) almost saturates the total string tension a JO 
The mechanism assumes condensation of the monopoles and, what is most important for 



1 



us, is borrowed from the electrodynamics of superconductors, i.e. is Abelian in nature. 
As a result, the dual formulation of the gluodynamics (see, e.g., || [|, ^ [| (jj) effectively 
reduces to the same U(l) formalism. The procedure is in two steps. First, one assumes 
Abelian dominance in the infrared limit of the gluodynamics itself. Then one constructs 
the effective U(l) theory a la Zwanziger keeping only these Abelian degrees of freedom. 
Thus, the crucial step is the Abelian dominance. There are various ways to introduce 
this dominance in the field theoretical language. In particular, the approach of Ref. [|J 
is to start with an effective theory for the dual gluons and to assume that it is similar 
to the Georgi-Glashow model. Thus, one postulates the existence of dual gluons, Higgs 
fields in the adjoint representation and the Higgs mechanism which reduces the symmetry 
of the dual worldQ from SU{2) to U(l). As a result, only the corresponding U(l) gauge 
field interacts with the quarks in the infrared. In other approaches [f| ^, |7| one assumes 
that in a particular gauge the masses of the ordinary gluons charged with respect to the 
remaining U(l) are larger than the mass of the "neutral" gluon. This "neutral" particle 
interacts both with color and magnetic charges as in the Ginzburg-Landau theory. 

The dual Ginzburg-Landau theories based on the Abelian dominance in gluodynamics 
are successful numerically as far as interaction of the heavy quarks in the fundamental 
representation is concerned. However, as was emphasized in Ref. ||13|| , the data on inter- 



action of the quarks in higher representations pose a very serious challenge to this picture. 
The point is that at intermediate distances (which cover in fact almost all the distances 
measured so far) the potential between higher representation sources consists of linear 
and Coulomb-like pieces. Moreover, the coefficients in front of the linear and Coulomb 



terms are proportional to each other fL4 | 



V T (R) « _T(T + 1)— ^- + T(T+l)aR, (3) 

where T labels the representation of the sources and a is independent of T. Then the 
crucial observation is that in an Abelian-dominant vacuum the quarks which are neutral 
with respect to the U(l) in question escape any interaction and there can be no linear 
potential like (|3]). This is true independent of the form of the interaction in the Abelian 
sector so that the Casimir scaling (|3]) seems fatal for the dual Ginzburg-Landau approach 
to the gluodynamics. More generally, Eq. (|3]) provides a crucial test for any model of 
confinement ||T3| . 



In this paper, we revisit the problem of constructing the dual gluodynamics by unifying 
the consideration of magnetic monopoles with \Q m \ = 1 and \Q m \ = 2. The point is that 
the fundamental monopoles with \Q m \ = 1 are introduced via the 't Hooft loop and 
are point-like in the continuum limit. To describe their interaction we need, therefore, 
understanding of the dual gluons on the fundamental level. The basic observation is 



that classically the fundamental \Q m \ = 1 monopoles correspond to U(l) solutions [fL6 
Therefore we introduce a dual gluon which is Abelian on the fundamental level. There is 
no assumption on the Abelian dominance whatsoever at short distances and the Abelian 



1 For simplicity we will consider SU{2) case throughout the paper. Generalizations to SU{N) are 
straightforward. 
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nature of the dual gluon is imposed on us by the well-known fact [[| [L6| that the monopole 
solutions of the Yang-Mills theory are the same U(l) solutions whose embedding into 
the SU(2) is arbitrary. We derive a Zwanziger-type Lagrangian which unifies standard 
gluons interacting with the color and a U(l) dual gluon interacting with external magnetic 
charges. Then it is a kind of a gauge freedom how to embed the magnetic U(l) into 
SU(2). There is no violation of the color since one averages over all possible embeddings 
or restricts oneself to evaluation of gauge-invariant quantities, i.e. the Wilson loops. 

We check the validity of this approach by evaluating the running of the coupling 
governing the Coulomb-like interaction of the external monopole-antimonopole pair. As 
expected, it does run as l/g 2 (r) , i.e. preserving the Dirac quantization condition. How- 
ever this is true only if the ultraviolet cut-off is introduced in a way consistent with the 
Dirac veto. We first postulate such rules and then check that they do follow directly from 
an analysis of the continuum analog of the 't Hooft loop operator worked out in [|T/J . 

Turning then to physics of larger distances, or physics of confinement we add to the 
Zwanziger Lagrangian a phenomenological piece ensuring Higgs mechanism in the dual 
sector. At this stage, one can assume that Abrikosov-like vortices with external quarks 
at the end points are formed. These solutions are pure Abelian of course but there is no 
violation of the color symmetry, the same as at the short distances. Whether the Casimir 
scaling (||) is observed or not becomes now a dynamical question since the string tension, 
&T is a function of the phenomenological parameters, that is of the masses of the dual 
gluon my and of the Higgs particle, m H : 

QJ 

V T (R) « -T(T+1)— ^ + T(T + l)<j T (m H ,m v )R. 

ix K 

We argue that the Casimir scaling holds in the London limit, m# 3> my. Thus, the 
approach developed does allow to reconcile the Ginzburg-Landau type mechanism of con- 
finement with the Casimir scaling. 

On the theoretical side, however, there is an important limitation on the use of the 
effective Lagrangian. Namely, upon accounting for the monopole condensation one cannot 
treat consistently quantum effects. The reason is that to impose the Dirac veto for virtual 
particles becomes now an intractable problem. Indeed, it was not simple already in case of 
radiative corrections to external monopole interactions (see above). And if the monopoles 
are condensed the Dirac strings are "everywhere". 

The paper is organized as follows. In Section we briefly review the classification of 



monopoles in non- Abelian gauge models 0, [16|, 18] . Following the construction of Ref . ]Tj 
in Section [5] we introduce the notion of \Q m \ = 1,2 monopoles in the continuum SU{2) 
gluodynamics. Then in Section |] the fundamental \Q m \ = 1 monopoles are considered. 
In particular, we obtain the interaction potential of the monopole-antimonopole pair and 
check the running of the coupling. In Sections |47^ we derive the Zwanziger-type 
Lagrangian for SU(2) gluodynamics and reconsider the \Q m \ = 1 monopole interaction in 
the new framework. Section || is devoted to discussion of \Q m \ = 2 monopole dynamics 
at large distances. We argue that within our approach the assumption of monopole 
condensation allows to reconcile the dual superconductor confinement mechanism with 
the Casimir scaling effect. Our conclusions are summarized in Section E. 
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2 Classification of the monopoles 



Since we are going to describe interaction of non-Abelian monopoles, the natural starting 
point is their classification. In the literature one can find a few classifications which look 
different. In this section we review two particular approaches, namely the dynamical 



classification of Ref. [[U| and topological one formulated in Ref. |yj (for review see ||). 
In the former case the monopoles are classified with respect to U(l) subgroups, while the 
latter approach utilizes the global structure of the gauge group, namely the non-triviality 
of the gauge group center. Therefore, we refer to them as U(l) and Z 2 classifications 
respectively. We conclude that the Z 2 classification applies in fact at short distances and, 
in our case, is relevant to the monopoles introduced via the 't Hooft loop. On the other 
hand, the U(l) classification is relevant to monopoles existing in the vacuum state of the 
gluodynamics. 



The Z 2 classification [|I^] is based entirely on topological arguments, according to which 
independent types of monopoles are in one-to-one correspondence with the first homotopy 
group of the gauge group. Therefore, there are apparently no magnetic monopoles at 
all in the SU(2) case. Nevertheless, for the gauge theory without matter fields in the 
fundamental representation (quarks) the actual gauge group is SO (3) = SU{2)/Z 2 and 
there exists a single topologically non-trivial monopole solution. We take the magnetic 
charge Q m of such monopole as unity. Note, however, that the charges Q m = ±1 are 
indistinguishable in fact. As for the charges \Q m \ = 2 they are equivalent, from this point 
of view, to no magnetic charge at all. The Z 2 classification might well be illustrated by 



Wu-Yang monopole solution JT9[. Indeed, the configuration 



A% = 0, At = 9^ £ akn^_ (4) 

g r z 

for Q m = 1 represents the (singular) Z 2 magnetic monopole , while for Q m = 2 it is gauge 
equivalent to the vacuum A = 0. 

Within the dynamical, or U(l) classification one looks for monopole-like solutions 
of the classical Yang-Mills equations. Where by the "monopole-like" solutions one un- 
derstands potentials which fall off as 1/r at large r. The basic finding is that there are 
no specific non-Abelian solutions and all the monopoles are, in fact, gauge copies of pure 
Abelian monopoles embedded into the SU(2) group. Moreover, using the gauge invariance 
one can always choose the corresponding U(l) group as, say, the rotation group around 
the third direction in the color space. According to this classification, the monopoles are 
characterized by their charge with respect to the U(l) subgroup and may have, therefore, 
charges, 

\Q m \ = 0, 1, 2, ... (5) 

Note that the topological solution constructed within Z 2 classification exactly corresponds 
to the \Q m \ = 1 monopole in the U(l) approach. For instance, the Q m = 1 configuration 
((D is gauge equivalent to 

a; = - <r* - (1 " CQS ^ , a* = a$ = o. (6) 

^ g 2r smt) 



4 



Therefore we can restrict ourself to the U(l) classification scheme only and consider 
every monopole configuration as a pure Abelian one embedded into the SU(2) group. The 
charge quantization condition (||) is then naturally interpreted as a classification of various 
admissible Dirac strings whose end points represent monopoles. Indeed, there is the usual 
requirement that the Aharonov-Bohm effect on the Dirac string is not observable. It 
turns out that the Dirac string corresponding to the \Q m \ = 1 monopole would be "seen" 
through the scattering of particles in the fundamental representation but is invisible to 
the gluons since their charges are twice as large. Contrary to that the \Q m \ = 2 string 
singularities cannot be detected by any test particle. It is worth emphasizing that in the 
standard lattice formulation of the theory the action is constructed in terms of infinitesimal 
Wilson loops in the fundamental representation, which implies in turn that the energy of 
the \Q m \ = 1 Dirac string is infinite in the continuum limit. Therefore, such Dirac strings 
and the corresponding \Q m \ = 1 monopoles can be introduced only as external objects 
via the 't Hooft loop. 

Within the Z2 classification, having the charge \Q m \ = 2 is equivalent to have no 
monopoles at all. This is not true, however, within the U(l) classification. As is empha- 
sized in Ref. ||, this apparent contradiction is resolved through the observation that the 
Abelian monopoles with \Q m \ = 2 are in fact unstable. The instability was demonstrated 
in Ref. || by analyzing the modes of the non- Abelian gluonic field in the monopole back- 
ground. Physically the instability corresponds to the gluon falling onto the monopole. 
Contrary to that the monopoles with \Q m \ = 1 have no negative modes and are stable in 
gluodynamics. 

Below we briefly describe the Q m = 2 monopole instability following Ref. 0. The 
quadratic fluctuation operator for the Q m = 2 magnetic monopole has a negative mode 
corresponding to the s-wave scattering of the charged gauge boson on the monopole. 
This negative mode corresponds to the instability for the charged boson to fall onto the 
monopole center due to attractive nature of the monopole-boson potential in the s-wave. 
The operator which determines the stability of the Abelian monopole solution is just the 
Schodinger equation for the a charged vector particle (gluon) in the field of the Abelian 
monopole ||: 

(a t . Sij - -^^ij) i>j = uj 2 4>i, (7) 

where ipi = ipi(r) are spatial components of the spin-one gluon wave function, uj 2 is the 
corresponding eigenvalue. The summation over the silent indices i, j, k — 1, 2, 3 is assumed 
here and below. The operator A r = — r _1 d 2 r is the radial part of the three dimensional 
Laplace operator and 

A ij = Q m S ij -r-(L-L - 5zk) , (8) 
where Sy = iEijkik is the gluon spin operator, L = ir x (d + 2iriQ m A) — Q m i/2 is the 
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angular momentum operator, A is the monopole field and r = r/r. The equations for the 
angular and radial parts of the wave function, ipi = Xi(9, <fi)R(r), are: 



A r - \\) R x (r) = uj 2 R x (r), (9) 



AijXjiB,^) = XXi(e,(f>). (10) 

For A < 1/4 there exists the continuous spectrum for Eq.(^), uj 2 > 0. This case 
corresponds to the scattering of the gluon on the monopole. In this case the Abelian 
monopole is stable in gluodynamics. According to Ref. || this situation is realized for 
the Q m = 1 monopoles. 

If A > 1/4 then there exists infinite number of the monopole-gluon bound states 
{oj 2 < 0), for the ground state u 2 = — oo. Physically this situation corresponds to 
the gluon falling down to the monopole and it is quite similar to the Callan-Rubakov 
catalysis |2IJ of the nucleon decay by monopoles. The monopole with charges Q m > 1 
always contain an eigenvalue A > 1/4 and thus they are unstable ||. 

In a much simpler way, one can demonstrate apparent irrelevance of the \Q m \ — 2 
monopoles by producing an explicit non-Abelian field configuration which looks as a 



\Qm\ — 2 monopole pair in its Abelian part but has no SU(2) action at all [T7j. Such a 



configuration is generated from the vacuum by the following gauge rotation matrix: 

Q = 



'ID 



where p is the azimuthal angle around the axis connecting the monopoles and Ad is the 
U(l) potential representing Abelian monopole pair: 



A^dx^ = 7:(—-—)d(p = A D (z,p)dp, (12) 



z± = z ± R/2 , p = x 2 + y 2 , r± = z± + p 2 . 

Note that the action associated with the Dirac string is considered in this case zero, in 
accordance with the lattice version of the theory (for details see ||17|| ). 

In this example the monopoles with \Q m \ = 2 are a kind of pure gauge field configura- 
tions carrying no action at all. Therefore, it seems natural to conclude that only the Z2 
classification is relevant. However, the crucial point is that \Q m \ = 2 monopoles do exist 
in the vacuum and, in particular, realize the equivalence of Q m = ±1 monopoles implied 
by Z 2 classification. Indeed, suppose that we start with, say, Q m = +1 configuration. 
Then a Dirac string carrying the Q m = —2 flux can be superimposed on this solution. It 
is important at this point that such a Dirac string costs no action (or energy). Then the 
radial magnetic field can also change its direction since it does not contradict the flux con- 
servation any longer. In a related language, one could say that the \Q m \ = 2 monopoles 
are condensed in the vacuum and that is why the magnetic charge can be changed freely 
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by two units. As a result, one expects that two \Q m \ = 1 monopoles would interact as a 
monopole-antimonopole pair, since their attraction corresponds to the lowest energy state 
of the system. 

Thus, any monopole-like configuration in gluodynamics may be considered within 
the U(l) classification scheme. At short distances, r <C Aq^ d only the fundamental 
\Q m \ = 1 monopoles can consistently be constructed. Moreover, they are infinitely heavy 
and thus only exist as external objects. On the other hand, \Q m \ = 2 monopoles do not 
produce any action at small distances, they are so to say empty at this scale. This fact 
implies in turn that the Q m = ±1 monopoles are indistinguishable from each other and 
therefore effectively describe the Z 2 monopoles. Contrary to that, the monopole dynamics 
at distances r > Aq^ d is much more subtle. In particular, one cannot prove any longer 
that the action associated with \Q m \ = 2 monopoles at large distances vanishes. We 
discuss the dynamics of \Q m \ = 2 monopoles in Section ||. 



3 Monopoles in the Continuum Limit 
of SU(2) Lattice Gauge Theory. 

We have shown that all monopoles in the SU(2) gauge theory may be considered in a 
unified way within the dynamical U(l) classification. It is worth emphasizing that the 
only possibility to have a consistent description of monopoles is to consider the continuum 
theory as a limiting case of lattice formulation. Indeed, we have argued that the classi- 
fication of monopoles reduces to the classification of all admissible Dirac strings, while 
in the conventional continuum theory the Dirac strings are not allowed at all. Thus we 
are encouraged to consider all admissible string-like singularities in the continuum limit 
of lattice gauge models. 

Following [17] we start from the SU(2) lattice gauge theory (LGT) with the standard 



Wilson action 

S = -J2S P , S p =- 2 TrU p , (13) 
p g 

where g is the bare coupling constant, p is an elementary plaquette and U p denotes the 
ordered product of link matrices in the fundamental representation along the boundary 
of p. The matrix U p = e lFp is defined in terms of the lattice field strength tensor F p = 
F p 2 a a /2. Therefore, the action (|J) depends only on cos(|F p |/2), \F P \ = (F^F*) 1 ' 2 and thus 
possesses not only the usual gauge symmetry, but allows also for the gauge transformations 
which shift the field strength by 47rfc, \F p \ — > \F p \ + Arck, k G Z. For \F p \ ^ we define 
n p = F p /\F p \ and let n a v be an arbitrary unit vector for \F p \ = 0. Then the symmetry 
inherent to the lattice formulation may be formulated as: 

F; - F; + Ann;, (n a p ) 2 = 1, e abc F h p n c v = 0, a = 1,2,3. (14) 

The symmetry ([14] ) which directly follows from the compactness of SU(2), is evidently 
lost in the naive continuum limit. Therefore, the action of SU(2) LGT in the continuum 
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limit is not given by the conventional expression ~ f{F® u ) 2 since it should respect the 
continuum analog of Eq. ([14]). 

One can show that in the continuum limit represents a singular two-dimensional 
string E® with a particular color orientation. The explicit construction of E" is as follows 
(see Ref. |17| for details). Let E^ denote a non self- intersecting surface with world-sheet 
coordinates x(a), parameterized by a a , a = 1,2: 

E^ = J d 2 a iiv 5 {A) (x - x(a)) , d 2 a^ = d 2 a h^ u , = e aP d a x^dpx v . (15) 

For given E M1/ consider the value of the dual field strength tensor *F° U = \e^ v \ p F^ p tangent 
to the surface: 

*F» = *F; V V . = d[»A° u] + g e^A^Al . (16) 

In close analogy with the lattice considerations above we define: 



n» = ^ for |*F| 2 = (*F a )V0, (17) 



while n a is an arbitrary unit vector for all points in which \*F\ = 0. Then the surface E" 
is constructed as follows: 

X% = J d 2 a llv n a 5^\x-x{a)), (18) 

where we have used the notations ([15]). With this definition the lattice symmetry relation 
( |T4P becomes in the continuum limit: 

F% - F; v + j*Z« u . (19) 

In order to construct a continuum action which respects Eq. (]19|) , note that this 
relation should be valid for any shape of E, Eq. (|1^). On the other hand, in the sector 
with no surfaces at all this action should reduce to the standard one 1/4 f{F® u ) 2 . The 
most straightforward approach would be to sum over all possible strings E: 

S(F) = -In J PSexp{ -\ j d 4 x (1£ + j *S^) 2 } . (20) 

Note that the Eq. (^) is much analogous to the Villain formulation of the compact elec- 
trodynamics. The corresponding partition function of the SU(2) gluodynamics becomes 
therefore: 

Z = J VAexp{ -S(F) } = J VAV^exp{ -i J d 4 x (J£, + j*Z%) 2 } . (21) 

To appreciate the meaning of the integral J "DE it is sufficient to consider one particular 
surface in Eq. (^1]). It turns out that E a exactly corresponds to the Dirac string associated 
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with \Q m \ — 2 monopole. Thus, the integration J PS implies invisibility of the \Q m \ = 2 
Dirac strings in the continuum limit of SU(2) LGT. The appearance of the additional 
color index is due to the Abelian nature of monopoles, which might have an arbitrary 
color orientation within SU(2). In particular, for non self- intersecting surfaces it is always 
possible to find a gauge in which all the Dirac strings are aligned along, say, third direction 
in the color space S a = 5 a,3 E. 

Note also that in Eqs. ( ^0| |2l| ) the integral j VE has, in a sense, a symbolic meaning: 
it is impossible to separate rigorously the measure £>E from the gauge degrees of freedom 
in T>A. A concrete example of this kind is provided by the Eqs. ( |TT| , |I~2"|) which show that 
the string singularities might be generated by gauge transformations. 

It is also possible to derive an analogous representation for \Q m \ = 1 monopoles. We 
have already noted that these monopoles are not inherent to the vacuum of gluodynamics 
and may only be introduced as an external objects via 't Hooft loop. Consider therefore 
the lattice definition of the 't Hooft loop operator: 

Huafa) = exp{-2 ^ S p } , (22) 

where Sj is an arbitrary surface spanned on a given contour j (the star symbol means the 
duality operation). In the path integral formulation the 't Hooft loop effectively changes 
the sign of the plaquette variables U p belonging to *£_,• : U p — > —U p , which can also be 
represented as: 

F p -> Fp + 2 ™ a p for p E , (23) 

where is defined exactly as before. The Eq. (|23|) is almost identical to Eq. ([14]) and 
therefore one can appreciate the continuum limit of fl2"2"|). Namely, the 't Hooft loop op- 
erator which introduces the \Q m \ = 1 monopoles into the vacuum of SU{2) gluodynamics 
is given by ]T7| : 

Hpj) = exp{S(F) - S(F (24) 

where the action S(F) is defined in Eq. (p0|). Note that if there was no explicit integration 
j DYj in ( p0| ) then the operator ( p4|) would insert the U(l) monopole pair in which the 
Q m = 1 and Q m = — 1 particles are different. This would be in clear contradiction with 
the general properties of the 't Hooft loop. Thus it is precisely the \Q m \ = 2 monopoles 
which ensure the Z 2 nature of the operator (|24]) . 

4 Z<i monopoles at short distances 

It is crucial for understanding the symmetries of the problem that the \Q m \ = 1 monopoles 
are point-like objects. Therefore, we are challenged to consider their interaction at ar- 
bitrary small distances or, in other words, at the fundamental level. Moreover, the 
monopoles in gluodynamics are in fact Abelian. Therefore, we immediately come to 
the paradoxical conclusion that the dual field, if any, is Abelian as well. There is no place 
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for a non-Abelian dual gluon because the monopoles do not constitute representations of 
the non-Abelian group. Then an immediate problem is how to match the non-Abelian 
symmetry of the gluon interactions with the color and the Abelian symmetry of the inter- 
action of the dual gluon with the external magnetic charges. We will discuss the problem 
in detail in this section. Basically, our approach is similar to that of Zwanziger JTJ] in case 
of electrodynamics. However, before going into this subject let us consider in more detail 



the expectation value of the operator fl24|) . 



Here, we will restrict ourselves to the rectangular T x R, T 3> R contours j. Then 

H(j) = (#(£,•)> ~ e~ TVm ^ R) (25) 

and we refer to V m m{R) as the inter-monopole (monopole-antimonopole) potential. For 
simplicity, we consider only the sector with no \Q m \ = 2 monopoles in which the action 
( pUP coincides with the standard one and therefore: 

H(j) = \\ VA ex P{ ~\ J d ^ + J * S ;V) 2 } • ( 26 ) 

It is easy to see that for the contours considered the potential V m fn(R) calculated with 
( p6|) coincides with the potential of a Z2 monopole pair. 

To derive V m m(R) at the classical level, note first that the gauge covariance of Eq. (|26|) 
allows to consider a particular gauge in which £!■ = S a,3 Ilj. Then the classical solution of 



the field equations corresponding to (|26|) is 



Aidx^ ee AyAxp, = — [ — dip, A 1 ' = , (27) 



1 






2^ 




r _ 




= x 2 


+ y 2 



z± = z± R/2 , p z = x z + y z , r| = z\ + p 2 



Moreover, it can be shown |T^, [TiJ that the general classical solution with minimal energy 
is a gauge rotation of (p7|) . Thus the classical inter-monopole potential is: 

V mrTl (R) (28) 
g 2 R 

It is worth emphasizing that the classical limit of the state created by the non-Abelian 
't Hooft loop operator is an Abelian (up to the gauge transformation) monopole-antimo- 
nopole pair. 

Next, consider radiative corrections to the Coulomb-like interaction ( ^8|) at short dis- 
tances. There are two specific problems which have to be addressed. First, the classical 
solution ([27|) has a Dirac-string singularity. Second, the radial magnetic field is also in- 
versely proportional to the coupling constant and at first sight this mixes up all orders of 
the perturbation theory. 

To sort the things out, note that the evaluation of radiative corrections addresses in 
fact two different, although closely related problems. These are the running of the coupling 
and stability of the classical solutions. Both aspects are unified, of course, into evaluation 
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of a single loop in the classical background. However, the running of the coupling can 
be clarified by keeping track of the ultraviolet logs, lnA^y alone and is universal since 
in the ultraviolet all the external fields can be neglected. Therefore, the coefficient in 
front of In Auv can be found by evaluating the loop graph with two external legs, i.e. the 
graph corresponding to the standard polarization operator in perturbation theory. This 
is true despite of the fact that the monopole field is of order 1/g. On the other hand the 
stability of the classical solution is decided by the physics in the infrared. Here one needs 
to consider the particular dynamical system, monopoles in our case, and the fact that the 
magnetic charge is of order 1/g can be crucial. 

Thus to derive the running of the coupling we have to consider only the terms quadratic 
in the external field and demonstrate that indeed the leading quantum corrections pro- 
duce g 2 —>■ g 2 (R). However, we have ignored so far the string singularity, which in fact 
invalidates direct perturbative calculations. Nevertheless, we argue that one can still use 
the perturbation theory supplemented by an additional rules, which are to remove the 
effect of Dirac strings so to say by hand. These additional rules and their meaning are 
discussed in detail in subsection 4.2. 

Consider the standard decomposition = + a M of the gauge potentials A^ into 
classical field A^, Eq. (|27D, and quantum fluctuations a M taken in the background gauge 
D^a^ = 0. Since the color orientation of is a functional of gauge potentials, in account 
of quantum perturbations the vector n a becomes 

n = n cl . + 8n, 5n a = e abc n b cL u c , (29) 

where the small parameter u is of order a M . But in fact one can neglect the deviation Sn 
and consider n = n c i instead (similar arguments may be found in Ref. ||). Indeed, the 
condition D^a^ = is invariant under a M — > a fl + D^X provided that the gauge parameter 
A satisfies D 2 \ = 0. The second order differential equation D 2 \ = supplemented by the 
boundary condition A = — ui on the surface always has a solution. Since n transforms 
in the adjoint representation the gauge transformation — > a M + -D M A brings the n a to 
its classical value n a cl = 5 a ' 3 . 

Then the effect of one-loop corrections to the potential (p8|) reduces to the calculation 

of 

Det 2 [ -V 2 ] ■ Det- 1 [ -V 2 5^ + 2igf flv ] , (30) 

where the first term is due to the ghosts, T>^ = <9 M — igA^ is a U(l) covariant derivative 
and ffj, u denotes the regular part of the classical field: 

F; v {A) = *Z jflu + f, v ). (31) 

In the leading log approximation the second term in Eq. ([30]) factories: 

Vet- l [-V 2 5^ + 2ig~f^} = Det" 4 [-P 2 ] ■ Def- 1 ^/^] . (32) 

Thus the paramagnetic interaction is free of string singularity. On the other hand, the 
diamagnetic interaction term apparently gives the contribution proportional to J (d^A^) 2 , 
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which includes also the string singularity. To decide whether the Dirac string becomes 
"visible" at the level of quantum corrections one can consider the evaluation of Det[— D 2 \ 



in the field of single monopole p2| , where the determinant may be calculated exactly. It 

f/J,u) 



turns out that Det[— V 2 } is proportional to the squared regular magnetic field, /(/^ 
and is string independent. 

Thus the perturbative and exact calculations of the leading InAj/y behavior lead to 
the distinct results, which differ however by the Dirac string contribution. The reason 
for this difference is not difficult to figure out. Indeed, in the theory of charged parti- 
cles interacting with monopoles there is a famous Dirac "veto" which forbids any direct 
interaction with string singularity. But exactly this prohibition is violated in ordinary 
perturbation theory, which uses the plane-wave basis in calculations. This is in sharp 
contrast with exact evaluation of Det[— T> 2 ], where one uses the basis of eigenfunctions 
which explicitly depend on given monopole background. The only way to make sense of 
perturbative result is to substitute d^A v ] — > thus discarding the string contribution. 
Evidently, this substitution goes beyond the field theory and should be considered as an 
additional rules in perturbative calculations. We will return to the discussion of this issue 
in section |4.2| . 

The analytical calculation of Det[— T> 2 \ in the case of monopole-antimonopole pair 



is impossible to perform f22|| . Nevertheless, in discussing the string independence the 
problem is in fact identical to the single monopole case. Thus the conclusion is that 
Det[— T> 2 \ depends only on the regular part of F£ U (A) from which the result for the inter- 
monopole potential at the one-loop level follows: 



In the next section we consider the expectation value ( |25| , |26D anew and introduce the 
dual representation for H(j). To test the validity of new formulation we calculate in 
section £l| the monopole-antimonopole potential up to the one-loop order and reproduce 
the running of the coupling fl33|). 



4.1 Lagrangian Approach 

In this section we derive a kind of dual representation for the expectation value (125 



Since our approach is close to that of Zwanziger M , it is worth to highlight the essentials 
of the latter construction. 

The Zwanziger Lagrangian which describes interaction of a U(l) gauge fields with 
Abelian point-like monopoles looks as: 

LzAAB) = l -{m-[dAA]) 2 + l(m • [d A B}) 2 + (34) 
+ l -{m-[dAA)){m-*[dAB\) - % -{m ■ [d A B]){m ■ *[d A A]) + ig e f-A + ig mJ m -B, 
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where j e (j m ) and g e (g m ) are electric (magnetic) currents and coupling, respectively, m M 
is a constant vector, m 2 = 1 and 

[AAB] ltt , = A ll B v -A v B ll , (m ■ [A A = m v [A A B]^ . (35) 

At first sight, there are two different vector fields, A, B introduced to describe interaction 
with electric and magnetic charges. If it were so, however, a wrong problem would have 
been solved because we need to have a single photon interacting both with electric and 
magnetic charges. And this is what is achieved by the construct (0). Indeed, the action 
( |34]) is not diagonal in the A, B fields and one can check that the form of the A, B inter- 
ference terms in fl3l}) is such that the field strength tensor constructed on the potentials 
A and B are in fact related to each other. To be more specific, let us consider A and B 
potentials which for example are due to a single magnetic current j m . Using the propa- 
gators ( |3~7|) to be derived shortly and the relation G + [m A (m • G)] + *[m A (m • *G)\ = 
valid for any antisymmetric tensor G^ one can show that 

[d A A] + i*[d A B] = -i(md)- x *[m A j m ] . (36) 

The right hand side of this equation is non-zero only on the Dirac string attached to the 
current j m and directed along m M . Thus the constraint (|36|) implies that there are only 
two physical degrees of freedom corresponding to the transverse photons which can be 
described either in terms of the potential A or B. Topological excitations, however, can 
be different in terms of A and B. 

The physical content of fl3H|) is revealed by the propagators of the fields A, B. Since 
the Lagrangian fl34[) is at most quadratic in A, B it is easy to integrate out potentials 
A, B and obtain the corresponding propagators. In the Feynman gauge they are 

(A„A V ) = (B„B„> = |f, (37) 

The propagators should reproduce, as usual, the classical solutions. And indeed, the 
(AA), (BB) propagators describe the Coulomb-like interaction of two charges and mag- 
netic monopoles, respectively. While the (AB) propagator reproduces interaction of the 
magnetic field of a monopole with a moving electric charge. 



Consider now the expectation value (25,26). Since it depends only on the current j 



and not on the particular shape of Sj one can take 

= J^lmM], (38) 

which is a particular solution to <9£j = j. For rectangular contours j which we consider 
the surface ( |3"5| ) has no self-intersection points. Therefore one can write 

£• = n a -^—[mAj) (39) 
J {mo) 
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thus extending the definition of the field n a into the entire space-time. Note that only the 
value of n a on the string Sj does matter while a particular way of extending n a (a) — > n a (x) 
is irrelevant. Then the expectation value H(j) is represented as: 

H(j) = \ jvA jvn exp{ -± j d 4 x (F a + jn a 7^> A j]) 2 } , (40) 

n 2 =l 

where the path integral is to be performed with a constraint 

£ abc n ^ F c, 1 *[ mAj ]) = o, a = 1,2, 3, (41) 
(mo) 

which is a consequence of Eq. (|T7[) . The constraint equation (II) may be implemented by 
an additional field x°"- 

H{j) = 1 j VnVAVx exp{ -1 1 d 4 x (F a + ^ n a 7^y> A j]) 2 + i ^-j 3 C}, (42) 



n 2 =l 



C M = C,( X ,n,F) = — £ * X «nV.*f) r (43) 



1 

(m<9) 

It is worth emphasizing that the field n a (x) is a kind of a fake variable in the representation 
fl42|). The path integral is clearly independent on n a away from the string (p8|), but at 
the same time n a (x) for x e £j is determined through flU]). 

The last step in deriving the dual formulation of fl25p6|) is to use the equality: 



exv {-\Jd 4 x(G + -^*[mAj]) 2 } = JvBexp{-Jd 4 xL(G,B,j)}, (44) 
L(G,B,j) = - A G 2 + l[ m .(dAB - i*G)] 2 + tjB, 



which is valid for any antisymmetric tensor field G^. The relation (m ■ G) 2 + (m ■ *G) 2 



G 2 /2 has been used to derive (fP]). Using Eq. (Kip and redefining B — > Z? — C we can 
finally transform (flOD into the form: 



H(j) = ^ VAVBVx Jvnexp{- J L(A,B, X ,n) + i— J j B } , (45) 



n 2 =l 



where Lagrangian L depends on the gauge potentials A a , Abelian vector field B^ and is 
independent on the current j: 



L = - A {F* V ) 2 + -[(m.[SAB- iW - e afec x V*F c ])J (46) 
Let us discuss the properties of the formulation (|45l) - (|4~6|) . 
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In absence of the external magnetic currents j the Lagrangian (pjf) reduces to the 
standard Lagrangian of SU(2) gluodynamics, as it should be. Indeed, for vanishing j the 
integration over B field cancels the last term in Eq. (f46|) regardless of n a (cf. Eq. (44)). 



The Lagrangian ( |4q ) as well as the interaction terms in Eq. ( f45|) obviously possesses 
SU(2) x U{1) gauge invariance: 



SU(2) : raV a = n -> fi-^nfi , U(l) : 5 M — > B fJi + d^a . (47) 



As is evident from the transformation law ([47]) the expectation value (|35|) may also be 
written as 

H(j) = JvnH n (j), (48) 



where H n (j) is given analogously to Eq. (|45|) but with a fixed n a . Thus -f^(j) is a gauge 
singlet part of an apparently gauge variant expression H n (j). Moreover, the integration 
n a is in fact redundant: the operator H n (j) is well defined for any choice of n a and 
coincides with H(j). This property of the 't Hooft loop is well known ||. In particular, 
in the Hamiltonian formulation the 't Hooft loop operator may also be defined in a gauge 
non-invariant way |23] . The matrix elements of this operator taken between physical and 



unphysical states vanish and this provides a well defined meaning to the 't Hooft loop 
despite of gauge variance. The same effect may clearly be seen in the formulation (|48|) . 
Indeed, since various choices of n a are related by gauge transformations one can show that 
Hm{j) = H n2 (j) = H(j). Note also that for every particular n a one still has to enforce 
the constraint (^T|), which for n a = 5 a ' 3 becomes: 

( Fa -7^*l mA j] ) = 0. a = l,2. (49) 
(mo) 

It is evident that the origin of the vector n a goes back to the freedom of choosing 
the color orientation of the monopoles. As is emphasized above the monopole solutions 
are Abelian in nature which means, in particular, that they can be gauge rotated to any 
direction in the color space. Thus, picking up a particular n a is nothing else but choosing 
a U(l) subgroup to which the monopoles belong. Note that the embedding of relevant 
£7(1), which is described by n a (x), may also be space-time dependent. In particular, the 
requirement (^) means that all the monopoles are in the diagonal subgroup of SU(2). 
We can either average over the directions of n a or fix it arbitrarily but evaluate only gauge 
invariant quantities. 

As far as the quantization is concerned, the Lagrangian (|46| ) reproduces the correct 
degrees of freedom of the free gluons. Indeed, in the limit g — > and for n a = <5 a ' 3 the 
Lagrangian fl46|) reduces to 

L° = ^[dAA a ] 2 + ^(m-[dAB - i*d A A 3 ]) 2 , (50) 
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which essentially coincides with that of Zwanziger (|34f) plus free Lagrangian for off- 
diagonal gluons. Quantization at this point is the same as in the case of a single photon. 
In particular, to the lowest order in g 2 the monopole-antimonopole interaction is given by 
a single 5-field exchange. Using Eq. ( p7|) one immediately reproduces the result (|28|): 

1 ?ir f Air 2 
H{j) = ex P {--(-) 2 y J,A-%} ~ exp{ T— A^(R) } (51) 

V mm {R) = -1.1 (52) 
^ R 

An apparent application of ([K]) would be evaluating the running of the coupling g 
in the expression fl5"2"|). And, indeed, exploiting the Lagrangian fl46|) one can approach 
this problem in a way quite similar to the case of pure electrodynamics, for a review and 
further references see ||. We consider the perturbative corrections to the potential (|5^) 
in the next section. 



4.2 Radiative corrections. 

In this section we consider the radiative corrections to the potential ([52]) starting from (f45|) 
with n a = <5 a ' 3 . Since for a constant n a the problem is quite similar to the consideration 
of Dirac monopole in the electrodynamics, there is not much specific about the derivation 
of the running of the coupling. And, indeed, our considerations overlap to some extent 
with those given in the original papers [22], p3L EH] and in the reviews O, |3| . Still, we feel 



that it is useful to present the arguments, may be in a new sequence, to emphasize the 
points crucial for our purposes. 

Therefore we first concentrate on the Dirac monopole with the minimal magnetic 



charge interacting with electrons and in one-loop approximation [g, |22], EJ], |25| . The 
crucial point which has been already mentioned above is that only loops with insertion 
of two external (i.e., monopole) fields can be considered despite of the fact that there is 
no perturbative expansion at all. Indeed, considering more insertions makes the graphs 
infrared sensitive, with no possibility for lnAjyy to emerge. 

Then, the evaluation of, say, first radiative correction to the propagator (B^B V ) in the 
Zwanziger formalism ([37D seems very straightforward and reduces to taking a product of 



two (AB) propagators and inserting in between the standard polarization operator of two 



electromagnetic currents. The result is [54 



(B^B u ){k) = 5j ^{l-L) + — 1— (^-m M m,)A (53) 
L = ^-\nA 2 uv /k 2 , 

where we have neglected the electron mass so that the infrared cut-off is provided, in the 
logarithmic approximation, by the momentum k. 

At first sight, there is nothing disturbing about the result fl53|). Indeed, we have 
a renormalization of the original propagator which is to be absorbed into the running 
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coupling, and a new structure with the factor (k ■ m)~ 2 in front which is non- vanishing, 
however, only on the Dirac string. The latter term would correspond to the renormaliza- 
tion of the Dirac-string self-energy which we do not follow in any case since it is included 
into self-energy of the external monopoles. What is, actually, disturbing is that according 
( |53D the magnetic coupling would run exactly the same as the electric charge, 

(A,A v ){k) = (54) 

violating the Dirac quantization condition. 

The origin of the trouble is not difficult to figure out. Indeed, using the propagator 
(AB) while evaluating the radiative corrections is equivalent, of course, to using the full 
potential corresponding to the Dirac monopole A^. Then, switching on the interaction 
with electrons would bring terms like Af, "fijif? . Since A^ includes the potential of the 
string electrons do interact with the Dirac sheet and we are violating the Dirac veto which 
forbids any direct interaction with the string. 

Let us demonstrate quantitatively that, indeed, it is the interaction with the Dirac 
string that changes the sign of the radiative correction. This can be done in fact in an 
amusingly simple way. First, let us note that it is much simpler to remove the string 
if one works in terms of the field strength tensor, not the potential. Indeed, we have 
H = H ra d + Hstring while in terms of the potential A^ any separation of the string would 
be ambiguous (see, e.g. Eq. (|27|)). 

Thus, we start with relating the potential, or energy to the interference term in the 
H 2 field: 

V mm = [ Hi-H 2 d 3 r. (55) 



Now, it is not absolutely trivial, how we should understand the product Hi ■ H 2 . Consider 
for example the magnetic field of a single monopole: 

H = H rad + H string , div H = , (56) 
1 1 

Hrad = ~ — grad (- r) , div H string = 5(r - n) . 

47r |r — ri| 

Then, by the analogy with the case of two electric charges, we would like to have the 
following expression for the interaction energy: 

V m m = / Hi ra dH 2 w d r = — — | r . (57) 

J 4tt |n - r 2 | 

Note, however, that if we substitute the sum of the radial and string fields for H 12 , then 
we would have an extra term in the interaction energy: 

f 3 1 1 
Vmm — / (Hi ra dH2 ,string + Hi string^-2 ,rad) d T = +2 - — . (58) 

J 4vr |n - r 2 | 

In other words, the account of the string field would flip the sign of the interaction 
energy! This contribution, although looks absolutely finite, is of course a manifestation 
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of the singular nature of the magnetic field H s t r ing ■ Note that the integral in fl58|) does 
not depend on the shape of the string. 

Thus, in the zero, or classical approximation we should write: 

Hi • H 2 = Hi jTOt 2 • H 2jrarf • (59) 

However, in the above example the electrons interact with full potential and therefore 
the first radiative correction would bring the product of the total Hi ■ H 2 which includes 
also the string contribution^. Indeed, the result in the log approximation would be as 
follows: 

,string + Hl,rad) ■ (H 2 

,string + H 2jra rf) — — LRipad ■ H 2)TO <z , (60) 



where at the last step we have used the observation (pq). 

Now, it is clear how we could ameliorate the situation. Namely, to keep the Dirac 
string unphysical we should remove the string field from the expression (|60| ) which arises 
automatically if we use the propagators (|37|) following from the Zwanziger Lagrangian. 
Thus, one should change H x H 2 in the expression ( pO"D to H l rarf -H 2 rad so to say by hand. 
The justification is that we should remove the effect of the string field from any observable. 
Therefore the result fl53f ) is inconsistent due to the interaction of virtual particles with 
Dirac string. The correct propagator [B^B V ) at the one-loop level is 

(B,B v ){k) = ^(l + L). (61) 

In particular, there is no Dirac string self-energy renormalization. 

One might wonder, how it happens that the couplings in the electric and magnetic 
potential run in opposite ways. Indeed, now we reduced the product Hi ■ H 2 to exactly 
the same form as the product Ei ■ E 2 in case of two electric charges (since the radial 
magnetic and electric fields are the same, up to a change of the overall constants). The 
resolution of the paradox is that the renormalizations of the electric and magnetic fields 
are indeed similar in the language of the Lagrangian. However, the small corrections to 
the Lagrangian and Hamiltonian are related as: 

5L = -5H. (62) 

Since E 2 and H 2 enter with the same sign into the expression for the Hamiltonian and 
with the opposite signs into the Lagrangian, Eq. ( p2] ) implies that the running of the 
couplings in the electric V e and magnetic V m potentials is opposite in sign. Which is, of 
course, in full agreement with expectations since V e ~ g 2 and V m ~ g~ 2 . 



Return now to the evaluation of first radiative corrections to the potential (|52|). For the 
sake of definiteness, we take n a = <5 a ' 3 and consider the Feynman gauge. The consistency 
of this gauge with Eq. (f49"D can be established in the same way as for the background 



2 At this point we assume in fact that Auy is smaller than the inverse size of the string, which is 
convenient for our purposes. Other limiting procedures could be considered as well, however. 
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a) b) 
Figure 1: The diagrams relevant for the calculation of one-loop correction to {B^B V ) propagator. 
The solid and dashed lines represent WW and ghost particles respectively. 



gauge condition above (see the discussion after Eq. (p9|)). Since Eq. (pE9|) is satisfied, the 
dependence on the field x a can be ignored^. Then the free field Lagrangian corresponding 
to (H) is: 

L° = -(m-[dAA}) 2 + -{m-[dAB]) 2 - i (m • [d A B]){m ■ *[d A A]) + (63) 

1 - 27T 

+ W^-d 2 )W^ + -[(OB) 2 + [dAf] + %cd 2 c + t(d 2 ( + trjd 2 r] + t — jB, 

where A = A 3 , W = 2- 1 / 2 (A 1 + iA 2 ), W = 2~ l / 2 (A l - iA 2 ) and c, c, C, C, V, V are 
complex-valued ghost fields. Furthermore, the interaction is given by: 

L in t = igA^ [W n (dW)-W^dW) + W u d ll W u -W v d ll W u \ + (64) 

+ 2ig[dAB]*Im(...) + glm(...) { [m A (m ■ [ d A A + i*d A B })} - [ d A A + i*d A B ] } + 

+ g [C(Ad)( - fj(Ad)rj\ + ig [c(Wd)C + c{Wd)r) - fj(Wd)c - C(Wd)c] + ... 

where Im(...) M „ = In^W^Wy) and dots denote terms quartic in the fields. 

The first quantum correction to the potential fl52| ) is due to the renormalization of 
(BB) propagator. As is evident from the interaction Lagrangian ( S4"D in 0(g 2 ) order 



only WW vacuum polarization effects should be considered. Therefore, we can use the 
Eq. ( p6|) to argue that the second term in the middle line of Eq. (|64] ) describes in fact the 
interaction of virtual WW pairs with Dirac string. In order to maintain the unphysical 
nature of the string such an interaction should be excluded as we have already discussed 
in detail. Thus we simply discard this term in the Lagrangian 



Then the evaluation of the radiative correction to the (B^B^) propagator reduces to 
the calculation of three diagrams shown on the Fig. 1, since vacuum loop with both BWW 
and AWW vertices vanishes identically. In the leading log approximation these diagrams 
give 

Fig. la = 5^ k 2 -f-k- (-8) , (65) 



(4tt) 



^ lb = ST k 2 H (-|) , (66) 



3 Note that this approach might be too naive and dynamics of x a should be considered more carefully. 
We leave this question for future investigation. 
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where we have neglected all the terms proportional to k^, k v and L = ln(A^y/fc 2 ). The 
(BB) propagator with account for the quantum correction fl66|) is therefore 

—i m a k T v g 2 L 2 2 i m a 'k T ' 

~ P ' 1 (4tt) 2 3 J + (4vr) 2 3 ' (kmy [ ^ ^ 

and is essentially of the form d53|) . As we have argued above the one-loop result of the 
type (0,0) is in fact wrong since it includes the interaction of virtual particles with 
string singularity. The correct treatment of Dirac sheet leads to Eq. (|6l[) instead of (|53|) . 
Therefore, the (BB) propagator at the one-loop level is 

Thus we have reproduced the running of the coupling in the inter-monopole potential 

The general conclusion which can be drawn from consideration of the radiative cor- 
rections is that the Lagrangian approach produces reasonable physical results only if it 
is supplemented by particular rules of handling the singular magnetic field of the Dirac 
strings attached to the monopoles. The rule is that we should always subtract the contri- 
bution of the strings, even if it is not vanishing within the perturbation theory. Essentially, 
the rule is that the product Hi ■ H 2 should not include the piece H ra d • Hstring- 

To justify this rule we have to go actually beyond the Zwanziger Lagrangian and 
consider a regularization procedure which takes into account the Dirac strings. It is 
worth emphasizing that within the lattice regularization the term H ra d • ^string is indeed 
absent. Let us recall the reader that on the lattice the Dirac string corresponding to 
\Q m \ = 1 monopole pierces negative plaquettes. This is true in the limit g 2 — > 0, or 
a — > 0. If one looks for small deviation from U p = —1, then the action does not contain 
terms, which are linear in perturbations. Note that this will be not true for the expansion 
around arbitrary U p ^ ±1. This result means in turn that in the continuum limit there 
is no term H ra d • Hstring- Thus, our naive removal of the effects of the Dirac string from 
the radiative corrections, is justified by the lattice regularization. 

So far we considered only graphs with two insertion of the external monopole field. 
At this level, the value of the magnetic charge is not important. For example, the same 
arguments would go through without change if we started with, say, monopoles with 
Q m = 2. But such monopoles are unstable || and this is a much more drastic effect 
than the would-be running of the coupling. The instability is revealed by looking for 
the negative modes in the external field, i.e. considering infinitely many insertions of 
the external field. In particular, the instability of \Q m \ = 2 monopole is due to the 
strong interaction of the spin of charged gluons with the magnetic field of the monopole. 
As a result of such interaction, the gluons fall onto the center of \Q m \ = 2 monopole 
thus invalidating any classical consideration. For the monopoles with \Q m \ = 1 there 
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(67) 



is no instability. Nevertheless, one could question whether the system of a monopole 
and antimonopole is also stable. There are no analytical ways to study the stability in 



this case. Numerical simulations |T7j show, however, no sign of instability, as one would 
expect. 

5 \Q m \ — 2 Monopoles and Effective Theories. 

So far we discussed the fundamental monopoles \Q m \ = 1 which can be visualized as 
classical infinitely heavy objects. Because of the infinite mass, these monopoles can be 
used only as probes of the QCD vacuum but play no dynamical role by themselves. The 
\Q m \ = 2 monopoles at small distances, as is mentioned above, are essentially trivial, 
since in the ultraviolet they are indistinguishable from singular gauge transformations. 
On the other hand, there exist very simple arguments that they can play dynamical role 
at distances r ~ ^qcd- Indeed, as far as the fundamental Lagrangian is concerned, 
the only role of the quantum corrections is the running of the non-Abelian coupling. 
In particular, if we consider a lattice coarse enough then gsu(2) becomes of order unity. 
Obviously, the same coupling governs the physics associated with any U(l) subgroup of 
the SU(2). However, if the coupling gu(i) becomes of order unity, then there is a phase 
transition associated with the monopole condensation pB| . Thus, it is very natural to 
assume that the monopole condensation occurs also in QCD since the running of the 
coupling allows to scan the physics at all the values of gsu(2) until one runs into a phase 
transition. 

The | Q m | = 2 monopoles condense in the vacuum and their physics is commonly 
described by field theories incorporating the Higgs mechanism [|], [5|, [| |7|]. This kind 
of models, however, fail to reproduce the Casimir scaling [^] and in this section we 
confront with this test our model with U(l) dual gluon matching color gluons in the 
adjoint representation 

In order to have a more quantitative description of the problem let us consider the 
partition function (fHf ) which by assumption correctly describes the SU(2) lattice gauge 
theory in the continuum limit. What is of particular importance is that in this represen- 
tation there is no much difference between (f2~T| ) and the expectation value (^), since for 
every configuration of S, Eq. fl2~T| ) essentially coincides with (|45|). Since it is known |T7 | 



that (£1|) depends only on the boundaries of S, we can directly apply the results of the 
previous section and represent (^) in the form: 

Z = Jvj VAVBVx j Vn exp{ - J L(A, B, X ,n) + i— j J B } , (69) 

n 2 =l 

where L(A, B,x,n) is given by (|46|) and J T>j denotes a well known measure of integration 
over closed world-lines of \Q m \ = 2 monopoles. 

The next step would be to perform the summation over monopole currents j. The 
integral j T>j may exactly be rewritten as a functional integral over complex scalar field 
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Jvjex V {-S(j)+ijjH}= jv^e W {~ + ^)$| 2 + V(|$|) ] } , (70) 

where is an arbitrary real- valued vector field. The concrete form of the potential V(|$|) 
depends crucially on the current-current interaction S(j). In particular, for S(j) = 
the potential vanishes, while more sophisticated current interactions lead to non-trivial 
V(|$|). For instance, the simplest form S(j) ~ j 2 results in the |$| 4 term in the effective 
Lagrangian for the field $. 

Unfortunately, the exact form of the current-current interaction and hence the poten- 
tial V(|$|) is so far unknown. Nevertheless, we may assume without loss of generality 
that in the infrared region the potential V(|$|) is of the Higgs type: 



V(\*\) = A(|$| 2 -77 2 ) : 



(71) 



where A and rj are phenomenological constants. Of course, the vacuum expectation value 
of the Higgs, or monopole field is of order Aqcd- Therefore, we get the following effective 
theory of SU (2) gluodynamics 



VAVBV& V X / T>n exp{ - / d 4 xL eS } 



(72) 
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+ -\(d + i— 5)$| 2 + K(|$|) 
* 9 



(73) 



which is presumably valid at low energies. 

Despite its apparent simplicity, Eq. (|73|) is highly speculative. Namely, it unifies point- 
like gluons described by the fields A a , B and $ which is an effective scalar field. One may 
justify the use of ( f73f) by assuming that the effective size of the monopoles with \Q m \ = 2 
is in fact numerically small, although generically it is of order Aqcd (see also p8|). 

What is also specific about the Lagrangian (|73|) is that the dual gluon is a U{\) gauge 
boson. The color symmetry is maintained by averaging over all possible embeddings of 



the (dual) U(l) into 577(2), see the discussion in section O . The confinement mechanism 
inherent to (|73| ) is the formation of the Abrikosov-Nielsen-Olesen string [2T3] which can 
be considered already on the classical level and which we discuss in detail below. More 
generally, the Lagrangian (|73|) exhibits the Abelian dominance in the confining region 
which is the dominance of Abelian-like field configurations in the full non- Abelian theory. 
This dominance is common to all the realizations of the dual-superconductor model of 
confinement and is strongly supported by the lattice data |ITJ. What we avoid, 
however, is the breaking of SU(2) to £7(1) which is inherent to the models j|] which start 
with the dual gluons in the adjoint representation and then add effective isospin one Higgs 
fields. 
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From the theoretical point of view, the most important limitation in the use of (|73|) is 
that it can be consistently treated on the classical level only. The difficulty to extend it 
to the quantum level is due to the Dirac strings. Indeed, perturbatively the Dirac veto is 
violated for virtual particles (see Section |4.2| ). When the monopoles condense the Dirac 
strings are filling the whole of the vacuum and there are no known ways to rectify the 
perturbation theory. In case of the Abelian Higgs model, it is even more convenient to 
use the dual language when the Dirac strings are attached to the electric charges. Then, 
to respect the Dirac veto one should impose the condition that the Higgs field vanishes 



along a line connecting charges [30]. The static Abrikosov-Nielsen-Olesen string satisfies 
this constraint [29|. However, there is no known complete set of solutions satisfying this 
boundary condition. 

To reiterate the point, in the standard field theory we can reproduce, for example, 
interaction of two static charges in two different but equivalent ways. Classically, we 
solve the equations of motions for the sources which are conserved currents. Quantum 
mechanically, we study the correlator (A^x) A u (y)) and substitute the vacuum zero- 
point fluctuations for the potentials A^ix). In case of the Abelian Higgs model with both 
electric and magnetic charges the quantum propagator cannot be constructed explicitly. 
Moreover, there is no doubts that it would involve the physics of strings. Thus, we are 
confined to the former approach when only classical static solutions are known. 

How far can we go within the classical approximation? Strictly speaking, the classical 
approximation by itself does not allow to make any physically interesting predictions. 
The point is that we are considering now the interaction of quarks which are degenerate 
in color. In particular, we are interested in the potential energy for heavy quarks in the 
color singlet state. This notion cannot be formulated entirely in terms of the classical field 
theory. Quantum mechanically, the projection to the color singlet state is implemented 
by considering the Wilson loop. It would be, therefore, desirable to formulate classical 
solutions in terms of propagators. 

Consider first the U(l) Zwanziger Lagrangian ( |34"D where the magnetic current is 



constructed on a scalar field which possesses also the Higgs-like self-interaction. Formally 



the (Afj, A v ) propagator can be easily found (see, e.g., |30| and references therein): 



1 Tfl 

(A^A V ) = — - — V—tf^ - m^) + ...), (74) 

At + rriy [kmy 

where dots denote the longitudinal terms. Taken at face value, the propagator ([74]) makes, 
however, no sense since it has a double pole in the (km) variable and thus depends explic- 
itly on the direction of Dirac string. All these pathological features manifest violation of 
the Dirac veto in the formal derivation of the propagator (|74|). The right way is to solve 
the classical equations anew supplementing them by the boundary condition of vanishing 
of the Higgs field along a line connecting the external charges ||30|| . However, one can 
introduce instead formal rules of integrating the double pole out and derive the potential 
energy in coordinate space corresponding to an exchange of a single A-field quantum 



In the London limit, uih ^> rriy the result coincides with the solution of the classical 
equations of motion. The tension of the string, which connects the external charges ±Ng 
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is: 

2 

7T HI 

a Ab (Ng) « iV 2 — m 2 y ln-f . (75) 

For an arbitrary relation between mjj and my one should allow for an adjusting overall 
factor ||. 

Although such a use of the propagator may look awkwardfj, it teaches us once again 
that the classical solution corresponds to a single quantum exchange. Note that the 
double exchange would result in infrared divergent expression. However, the rule is that 
one should not iterate the quantum exchanges to reproduce a classical solution (for a 
similar analysis of the instanton solution see p^). 

Next, consider classical solutions corresponding to the Lagrangian ([731). Clearly, it 
has the same, Abelian-type solutions for any choice of the vector n a . We can therefore 
supplement the propagator ( |74"D by the Kroneker symbol 8 a ' b in the color space and reduce 
the effect of the classical solutions to a single quantum exchange. Or, we can fix n a = 5 a ' 3 
and consider the Wilson loop in the single-quantum exchange approximation. Note that 
all these remarks apply also to a single-gluon exchange as well. In this sense, there is no 
much specific about the Abrikosov-Nielsen-Olesen string exchange. Thus the criticism of 
Ref. [|13| does not apply in our case. We do have a U(l) dual gluon, but no U{1) subgroup 
is singled out for the ordinary gluons. The fixation of U(l) is a matter of gauge fixing 
and convenience. 

What is specific for the one-gluon exchange is that the overall coefficient is fixed by 
the condition that the gluon is coupled via the generators of the gauge group. In case of 
the string, the string tension which plays the role of the coupling is a dynamical quantity 
depending on the parameters of the theory, that is m#y. In this way we derive linear 
potential for heavy quarks in any representation T. As for the coefficient in front of the 
linear rising term, it should be now calculated explicitly as function of my, run'- 



1 T 

— M2T 3 ), (76) 



(JT 

2T+ _ 

T 3 =-i 

where a^b denotes the string tension, calculated in pure Abelian Higgs model with the 
same values of w,h,v an d for external charges of magnitude ±2T 3 . In the London limit, 
vtlh ^> my, the string tension a"Ab(2T 3 ) ~ T 3 2 and we reproduce the Casimir scaling since 

T 

1 „o 1 



£ Tl = -T(T + 1) 



2T+1 ^ 6 3 

T 3 =-T 

So far we discussed only embedding of the U(l) into SU(2). As far as the physics of 
large distances is concerned, there is a subtle point how to define the U(l) dynamically. 

4 Alternatively, we could use the language of quantum mechanics. Indeed, our problem is similar to, 
say, the problem of hyperfine splitting. The role of the Wilson loop is then to ensure that the external qq 
pair is in the T = state. While the knowledge of the classical solution allows to evaluate the potential 
in case when the sources have a definite (and opposite) T3. Averaging of the potential over the wave 
function gives the final result. 
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As is mentioned above the £7(1) arises in lattice simulations through gauge fixing [11 



The monopole properties depend on the particular choice of U(l). The Maximal Abelian 
projection turns most successful as far as the interaction of the heavy quarks in the 



fundamental representation is concerned |10|. 



It is plausible to assume that the Lagrangian (|73|) with n a = 5 a ' 3 describes the data in 
this gauge, although at present there are no theoretical arguments in favor of this. More- 
over, if this assumption is valid, then a new phenomenological problem arises. Namely, 
the structure of the confining string for the quarks in the fundamental representation is 



best described in the Bogomolny limit |33| , itlh ~ ray. While the Casimir scaling can be 
reproduced only when m# ^> my. Thus, an overall fit can be obtained only at price of 
compromising the quality of the fits. Which might be not surprising since we are relying 
on the classical approximation to fit the data obtained via the lattice simulations which 
are fully quantum. 



6 Conclusions. 

The main objective of this paper was to develop an approach to gluodynamics which 
would allow to describe interaction both with color and magnetic charges. We considered 
first external magnetic charges which can be introduced via the 't Hooft loop. We used the 
Zwanziger-type approach which formally starts with two different fields for ordinary and 
dual gluons, interacting with the color and magnetic charges, respectively. The doubling 
of the degrees of freedom is avoided then by imposing constraints. Our central point is 
that it is natural to combine ordinary gluons in the adjoint representation with £7(1) dual 
gluon. The reason is that the external monopoles are in fact Abelian objects and do not 
constitute representations of the non-Abelian group. The fixation of the £7(1) is a kind 
of new gauge freedom. In particular it can be fixed to be space-time independent. The 
overall non-Abelian symmetry is maintained either by averaging over all embeddings or 
by confining the consideration only to gauge-invariant quantities, i.e., the Wilson loops. 

We checked the viability of the approach by considering the running of the coupling 
governing the interaction of two external monopoles at short distances. We have shown 
that the coupling runs as expected, i.e., the Dirac quantization condition is preserved with 
account of the radiative corrections. To derive the running of the coupling it is essential 
to respect for the virtual particles the Dirac veto which forbids any direct interaction with 
the Dirac strings attached to the external monopoles. In the field theoretical language 
the veto is manifested through a particular ultraviolet regularization. This regularization 
can either be postulated or derived by analyzing the continuum analog of the 't Hooft 
loop. 

We turned then to consideration of large distance physics which is commonly described 
in terms of monopole condensation. Here, the monopoles with a double charge are meant 
which are not point-like, unlike the external monopoles introduced via the 't Hooft loop. 
What unifies the description of the both types of the monopoles is that the dual gluon 
is a £7(1) boson in both cases. This does not imply any violation of the color SU(2) 
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either at short or large distances. Classical, Abrikosov-type solutions to the effective 
theory produce linear potential for heavy quarks in any representation of the color group 
and allow, in particular, to incorporate the Casimir scaling. The use of the effective 
theory is limited, however, by the necessity to impose the Dirac veto for virtual particles. 
In practice, there is no known way to go beyond the classical approximation. At the 
quantum level, the physics of strings seems to be indispensable. 
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